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Abstract 

A free massless scalar field is coupled to homogeneous and isotropic loop quantum 
cosmology. The coupled model is investigated in the vicinity of the classical singu- 
larity, where discreteness is essential and where the quantum model is non-singular, 
as well as in the regime of large volumes, where it displays the expected semiclassical 
features. The particular matter content (massless, free scalar) is chosen to illustrate 
how the discrete structure regulates pathological behavior caused by kinetic terms 
of matter Hamiltonians (which in standard quantum cosmology lead to wave func- 
tions with an infinite number of oscillations near the classical singularity). Due to 
this modification of the small volume behavior the dynamical initial conditions of 
loop quantum cosmology are seen to provide a meaningful generalization of DeWitt's 
initial condition. 



1 Introduction 



Since the early days of canonical quantum gravity |l|, isotropic cosmological models have 
been popular test objects. Due to the symmetries their number of degrees of freedom is 
finite and small so that they become accessible to quantum mechanical methods. In the 
geometrodynamical approach quantum states are considered to be continuous functions of 
one or more variables, typically metric or exterior curvature quantities of space-like slices, 
and of matter variables if the system is coupled to matter. 

In the classical limit at large volume such models are well behaved and reproduce their 
classical counterpart. The problem of the classical singularity at zero volume, however, 
is not cured, which is related to the fact that spectra of geometric operators remain con- 
tinuous in standard quantum cosmology. (Nevertheless, at the level of expectation values 
the expansion/contraction velocity, for example, may remain finite and discrete features 

*e-mail address: bojowald@gravity.psu.edu 
t e-mail address: franz@piLysics.muni.cz 



1 



may emerge M.) On the other hand, loop quantum cosmology M inherits from quantum 
geometry a discrete structure of geometry [||, [| which is most relevant at small scales. 
Still, at large volume standard quantum cosmology emerges as an approximation J7| which 
is very good in the semiclassical regime but not applicable near the classical singularity, 
where quantum effects of gravity become dominant and the discreteness of geometry at 
the scale of the Planck length becomes essential. The discreteness leads to a resolution 
of the singularity problem through the following mechanism: Already at the kinematical 
level, there is an indication for a natural curvature cut-off since the classically diverging 
quantity a -1 has a finite quantization with an upper bound of the size of the inverse Planck 
length |8j. Furthermore, an investigation of the dynamical evolution confirms that the clas- 
sical singularity does not appear as a boundary but instead allows a well-defined evolution 
through it ||. 

In the context of the present paper we are mostly interested in another consequence 
of loop quantum cosmology, namely that it predicts dynamical initial conditions which 
are derived from the evolution equation |1(J. The issue of initial conditions in quantum 
cosmology has been widely discussed [pi [II], [T2|, in part as an attempt to deal with the 
singularity. De Witt's initial condition, which is closely related to the outcome of the dy- 
namical initial conditions, requires the wave function to vanish at the classical singularity 
(which, however, with continuous geometrical spectra does not solve the singularity prob- 
lem since one can still get arbitrarily close to vanishing volume fljj]). On the other hand, 
it is well known that DeWitt's condition is not applicable in general because in most cases 
it would predict an identically vanishing wave function. A particularly thorny issue, also 
for other boundary proposals, is the fact that solutions to the Wheeler-DeWitt equation 
often exhibit an infinite number of oscillations between vanishing scale factor and any finite 
value. In such a case, the limit of the wave function for a — > is not always well-defined 
and one cannot even impose initial conditions at a = 0. The origin of the infinite number 
of oscillations is the kinetic term in a matter Hamiltonian which is proportional to a~ 3 . 
At small a, the Wheeler-DeWitt equation requires the derivative of the wave function to 
be proportional to the square root of the kinetic term, which diverges for a —>■ 0. Usually, 
one tries to avoid this problem by choosing a wave function which is independent of the 
matter field at small a so that the matter momentum vanishes and the kinetic term is 
identically zero. However, for larger volume the wave function must depend on the matter 
field non-trivially which forces a dependence at small volume, too. Even if the dependence 
is only weak, the diverging kinetic term will eventually dominate when a is small enough. 
In the present paper we investigate if there is a more natural way to deal with this problem 
from the point of view of loop quantum cosmology. 

Since this problem is caused by the kinetic term independently of the potential, we will 
analyze it in the most simple setting, which is isotropic loop quantum cosmology coupled 
to a free massless scalar field. In this way, the total quantum Hamiltonian acts on the state 
function containing one degree of freedom of geometry and one of matter. The formulation 
is based on |3J |8| where a discrete orthonormal basis of states of homogeneous and isotropic 
quantum geometry has been established which provides eigenstates of the volume operator. 

This paper is organized as follows. In section 2 the quantum Hamiltonian constraint 
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equation with a massless scalar field is formulated, leading to a difference equation which 
is solved in section 3. In section 4 the continuum and semiclassical limit are considered 
and possible physical consequences are discussed. Section 5 deals with a comparison of 
different proposals for initial conditions. 



2 The quantum Hamiltonian constraint 

We begin with the Lagrange density of a massless scalar field on a spacelike slice of a 
canonically decomposed 4-dimensional space time manifold, 



r — 1 

*-<f> 2 



2 " (V# 



where V denotes the spatial derivative on the slice and the dot the derivative w. r. to the 
time coordinate. In our case, the spatially flat Friedmann model with the metric 

ds 2 = -dt 2 + a 2 (t) [dx 2 + dy 2 + dz 2 ] , 

the total matter Lagrangian becomes 

d 3 i^£ = \ I d 3 xa 3 2 , 
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where, in order not to disturb homogeneity, is assumed to be spatially constant so that 
the gradient term vanishes identically. The integral over the coordinate depending term, 
i. e. the total volume divided by a 3 , will be set equal to one which can be achieved by an 
appropriate compactification. For the Hamiltonian of the field we obtain 

H^IpIcT 3 (1) 

where p^ = o?4> is conjugate to <fi. 

Together with the gravitational part of the Hamiltonian constraint we obtain in the 
isotropic case 

H = -Q 1 ~ 2 K- 1 c 2 ^\p~\ + = -Q^k-'c 2 ^ + \\p\~ipl = (2) 

where the canonical gravitational degrees of freedom are the isotropic connection compo- 
nent c and (densitized) triad component p which fulfill {c,p} = |7K. Here, k = 8ttG is 
the gravitational constant and 7 the Barbero-Immirzi parameter [|l(| whose value does not 
affect the classical behavior (but it is important for the quantum theory where it controls 
the continuum limit). Because p is a triad component which has two possible orientations, 
it can take both positive and negative values. This is also true for the scale factor a which 
is the isotropic co-triad component and related to p by p = sgn(a)a 2 . From now on, how- 
ever, it will be sufficient to consider only positive p and a (though we keep absolute value 
signs in some formulae for the sake of generality). The connection component c is related 
to the extrinsic curvature and therefore to the time derivative of a by c = \a. 
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In the absence of a potential, the matter momentum p^ = 00 is a constant in time and 
the solution 

c=^^/k/3uj/p (3) 

of (0) only depends on p. To understand the classical evolution in a coordinate time t 
(with lapse function N = 1) we compute 



p = {p, H} = A'j x c^fp = 2^k/3lup 2 

yielding p(t) = a(t) 2 = (y/3Ku(t — to))i, i.e. an eternally expanding universe. 

Now we quantize the field canonically by assuming a wave function x{4>) an d the canon- 
ical momentum acting derivative operator on it, 

h := (4) 
so that the Hamiltonian for spatially constant fields becomes 



d 
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H * = -W\*\ tyi ^ 

(it is denoted since we will later introduce another operator in which also a is 
quantized). 

For a quantization in the complete system of gravity and matter we need a quantization 
of the inverse power \a\~ 3 which diverges classically close to the singularity. In standard 
quantum cosmology this would simply be quantized to a multiplication operator acting 
on wave functions depending on a, which does not cure the divergence. Loop quantum 
cosmology, on the other hand, can easily deal with this problem: while the volume operator 
has zero eigenvalues and so no well-defined inverse, there are well-defined quantizations of 
inverse powers of the scale factor ||. This inverse power of the scale factor is essential for 
the issues studied in the present paper. The details of the quantization of the matter field 
is irrelevant; one can also use quantization techniques inspired from loop quantum gravity 
(see e.g. fOp). 

To specify the gravitational part of the wave function in loop quantum cosmology we 
start in the connection representation where an orthonormal basis is given by j|] 

(c\n) = r V2 / , neZ. (6) 
V2 sin |c 

These states are eigenstates of the volume operator 

V\n) = (| 7 ^) f V(M-1)M(M + 1) |n> == ^(M-i)N (7) 

and of the inverse scale factor operator (we only use the diagonal part of the operator rhjj 
derived in 0) 

[Sp|n> = 16( 7 ^)- 2 (^V^ - y/Vi^Yln) (8) 
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with lp = v nh being the Planck length (in (^j) VI i is understood to be zero). As discussed 



in [17]], the quantization of the inverse scale factor is affected by quantization ambiguities. 
The effect of different choices will be discussed later; they do not lead to substantial changes 
in most of the following calculations and results. 



The action of the gravitational Hamiltonian on the basis states \n) is P, [13 



H grav \n) = 3 7 - 2 («7£) 1 sgn(n) (v |w - V\ |nhl ) (\n + 4) - 2\n) + \n - 4» . (9) 

Now we write the states of the coupled system matter plus gravity as a superposition 
of geometry eigenstates 

oo 

\*)= E s ^)i™> ( 10 ) 

n=— oo 

with ^-dependent coefficients s n (<f>) which represent the state in the triad representation. 
We use (§) to quantize the inverse volume in the matter Hamiltonian (HI). Since the 
resulting operator is diagonal in the basis states \n) we can define 

H^n)® |0) =: |n)®i^(n)|0) (11) 

for each \n) and arbitrary \<p). For the massless field we only have to insert the eigenvalue 
of |a| _1 in a state \n), so we obtain 



H,(n) = -|n 2 16 3 ( 7 ^)- 6 ^V |w - ^ |Bhl J ^. (12) 

Finally, a state |s) is annihilated by the total Hamiltonian, the sum of H gr3iV in (P) and 
Hfj, in (|ll|) with ([T2|), if s n fulfills (we absorb the sign of n appearing in (|j) into the wave 
function) 

(Vi| n+4 | - V|| n+4hl J s n+4 (0) - 2 (Vi w - Vi| nhl J s n (4>) (13) 
i) s n _ 4 (0) = ah 2 - y / %^I) -jrs s n (0) 



+ ( Vl |n _ 4 | - Vl|„_ 4 | 



with 

This is a difference equation of order 8 for the wave function s n (<f>) in the internal time 
n. A priori, there could be a problem at zero volume since some of the coefficients of the 
difference equation can vanish due to V_i = V_i = Vq = 0. In fact, the coefficient of so{<p) 

- and only this one — always vanishes (note that, unlike in a classical equation, |a| _1 
annihilates the zero- volume state which represents the classical singularity). Thus, there 
is always a solution s n ((f>) = s o (0)5„o which is orthogonal to all other solutions and need 
not be taken further into account. The only free function then is s n (4>) for n = 4, from 



5 



which one can compute sg, and so on, yielding all coefficients S4 m (0) as functions of s 4 (0). 
The intermediate values S4 m+ j(0) for % = 1, 2, 3 are essentially fixed by the principal series 
S4m(0) by requiring that the wave function does not vary strongly on small scales (i.e. that 
it is pre-classical fllQl). 

The simplest choice for the free function 54(0), besides a constant, is an eigenfunction 
of the matter Hamiltonian -f^,(4) for n — 4, 

s 4 (0) = X (<f>) ■= e***/* (15) 

where u is the eigenvalue of p^. (A spatially constant massless field on Minkowski space 
with a non-vanishing energy eigenvalue would not be possible, but here geometry acts as a 
potential.) Because x(4>) is an eigenfunction for all H^n), the ^-dependence of all the s 4m 
is the same. This will simplify the analysis since we can compute the wave function from an 
ordinary difference equation rather than a difference-differential equation; the combination 
s n ((p)e~ luJ ^^ h will be ^-independent. 

We do not intend to justify the initial condition (|T^) with any given u by a physical 
argument; rather, we choose it to simplify the subsequent calculations while retaining 
qualitative aspects as to the behavior of the wave function s n (<f>) close to the classical 
singularity in the presence of matter. The advantage of an initial condition like flT5|) is 
that the gravitational and the matter degree of freedom separate such that the n-behavior 
of the wave function is still given by an ordinary difference equation. If a more complicated 
function than ( |l"5"|) is required, it can always be constructed as a suitable superposition of 
our solutions with different u. 



3 Solutions of the Hamiltonian constraint 

Having made the above observations about Sq and S4, we begin with inserting n = 4 into 
([jjD and obtain s$ in terms of s 4 = x(</>), 

W) = T7 77 X[4>) 

V4 — I/3 

and, in consequence, all 

/ / \ 2(V2m-2 — Vims) ~ a ^ 2 (V^2m-2 ~ V^m-s) 6 / ,s „m 

S4m(0) = 77 77 S 4m - 4 (0) (16) 

SAm-sW) tor mEN. 



2m ~ V2m-1 



For an explicit calculation of the coefficients s^ m it is convenient to introduce some 
shorthand notations. We define 

D n := V n - V n -t and W n := (^/V n - ^V^f , 
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then we split off the ^-dependence by defining ^-independent coefficients r m by r = 0, 
ri = 1 and 

SAm((f)) =: ■= — =r— ^ p- x(0) for m > 2 . 

Inserting this into (|16| ) yields a recurrence relation for the r m : 

(2D 2m -2 - au 2 W 2m - 2 ) — D 2m -2D 2m -i ■ r m _2- (17) 

An explicit solution can be found after calculating a few r m , and confirmed by induction: 

m—l 

r m = ^(-a^ 2 ) fc ^ j'i(j 2 - ji) • • ■ (jfe - ifc_i)(m - j k )W 2jl ■ ■ ■ W 2jk D 2jk+1 ■ ■ ■ D 2jm _ x . 

(18) 

This formula has the following meaning: j\,---,jk are integers between 1 and m — l, 

(Tfi 1 \ 

, j possibilities to choose k numbers out of {1,2, ... ,m — 1}, 

arranged in increasing order. The numbers jk+i, ■ ■ ■ ,j m -i appearing as labels of D n are 
given by the remaining values in {!,..., m — l}\{ii, • • • , jk} i n & n arbitrary ordering. 
The case k = is to be understood in the way that from the factors containing j there 
remains only the last one, m — j , with j defined to be zero, so that this contribution is 
mD 2 ■ ■•D 2 (m-i)- 

For the purpose of numerically calculating the coefficients in terms of a given initial 
one the recurrence relation ([17]) may be more convenient than the solution (|Tg|). 

Inserting n = into ([13]) leads to s_ 4 (0) = s 4 (0). Further on, (|I3|) is symmetric 
in the sense that the relations between s_ 4 , s_ 8 , s_ 12 ,. . . are the same as the respective 
relations for s 4 , sg, si 2 ,. . . and so the series s 4m is symmetric in m. By the argument of 
pre-classicality [0] this symmetry applies to the remaining series s 4m+ j as well. With the 
physical interpretation of n as an internal time we have obtained time symmetry of our 
wave function with respect to the classical singularity. 



4 Continuum and semiclassical limit 

For large n, when V n becomes very large in comparison with lp, the discreteness should only 
lead to small corrections, i. e. the discrete time evolution from n to n+4 should become well 
approximated by a continuous evolution. Under this assumption the difference equation 
should be approximated with high accuracy by a differential equation for a continuous 
wave function ip(p,(fi) = s n(j ,)(0), where n{p) = 6p/7/p can be derived from the volume 
eigenvalues (0). As discussed in [|10|| , among the solutions of the difference equation there is 
always one which is slowly varying at small scales for large n and so justifies the assumption 
of almost continuous behavior. 
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4.1 Large volume behavior 



To approximate ([13]) by a differential equation it is convenient first to make a slight change 



of variables. We define the new variable (working only with positive n > 0) 

t n (<f>) := (jll)- 1 (Vs - Vn-x) s„(0) (19) 

which for large n is approximated by t„(0) ~ ^^/p(n)tp(p(n), (f)) =: tjj(p(n), (f)). In terms of 
t n (4>) the evolution equation becomes 



(v^-v^) 6 d 2 



t n+4 (0) - 2t n (0) + t n _ 4 (0) = «r '— — t n {(j>). (20) 

l/n — V 2_ X Q</) 

Furthermore, we need the asymptotic approximations 

y y „ 24-*( 7 Z£)§7i* = i 7 Z!v'p(n) 



(^Vf-V^) 6 - (ife) i (7$)*n- i = (h^) 8 P(n)- f - 

For the differences operator on the left hand side of (PO) we obtain in the continuum limit 
16d 2 /dn 2 = |( 7 Zp) 2 9 2 /9p 2 , acting on the function ip(p,(p), so the asymptotic equation 
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becomes 

~ Alt , <9 2 ^ 



L P 2 

■v 



3nh 2 *" dp 2 d(f) 2 \^'^ ) '' ( ^ ^ 

Note that a standard quantization in geometrodynamics would have ordering ambiguities 
in the gravitational part of the Wheeler-DeWitt equation, whereas the derivation here 
leads to a unique ordering in loop quantum cosmology [Q. 

The variables in this partial differential equation can be separated by a product ansatz, 

#(p,<f>) = N(p)x(4>) , 

so that we obtain 

4/p 2 N " X" 2/ t 2 

■p —— = — = — uj jn = const. 



3Kh 2 ^ N X 

We assume uj 2 > 0, because the matter Hamiltonian is expected to have a positive spec- 
trum. With this assumption the solution \ coincides with the function used for s 4 in fll5|), 
denoted by the same letter. The ordinary differential equation for N(p) has the solutions 
N(p) = p x with 



A = |±Vl-!^V- (22) 

The asymptotic differential equation (pip applies also for t im+i for i — 1, 2, 3, so, in order 
to have a smooth function for large n we assume the 0-dependence of the remaining three 
free coefficients si, #2 and S3 to be the same as that of S4, given by the function x(0), 
endowed with constant factors which could be determined by a numerical analysis. 
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m 

Figure 1: The 4>- independent part S4 m ((j))e~ tul ^^ h oi the discrete wave function (x) compared 
to a continuous solution to the Wheeler-DeWitt equation (dashed line) for u 2 = 500a;;?- 
Strong deviations occur only at very small n = 4m (see Fig. |j) and will be discussed later. 

The last result fl2"2"|) deserves some discussion, because it provides a threshold for a 
qualitatively different behavior of the solution N(p), according to the sign of the expression 
under the square root. If u is smaller than a critical value, determined by 

u 2 = iz£/« = \nh 2 , (23) 

the exponent A is real and for uj going to zero one solution for N(p) approaches p, the 
asymptotic function for the vacuum calculated in [[|, the other one approaches a constant. 
For u larger than u c 

N(p) = p^p ±iQ = p l 2 e ±M**P (24) 

with 

n = l^3nu% 4 - 1 = i>r V3^7« - h 2 . (25) 

Reconstructing if)(p, (f>) from ip(p, <fi) we finally obtain 

tfj{p,<f)) = 2e ±imosp e iuj ^ /h . (26) 
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Formally this solution is a quantum mechanical wave function of a particle moving from 
small values of p to larger ones with decreasing speed, whereas solutions with real powers 
of p, coming from undercritical values of uj, are lacking such a dynamical interpretation. 
An example of an oscillating solution can be seen in Fig. |l] for uj 2 = 500^, a non-oscillating 
one with uj 2 = \uj 2 c in Fig. 0. 




5 10 15 20 25 30 35 40 

m 



Figure 2: A subcritical {uj 2 = \uj 2 ) wave function s im {(j))e~ lu "^^ h compared with the con- 
tinuous DeWitt wave function (solid line) and a continuous wave function which does not 
fulfill DeWitt's initial condition (dashed line; any such solution diverges at a = 0). 

Note that the presence of a critical frequency uj c depends on the ordering of the Wheeler- 
DeWitt equation; choosing, e.g., the ordering pd/dp{pd/dp) instead of p 2 d 2 /dp 2 would give 
oscillating solutions for all non-zero uj. Therefore, standard quantum cosmology cannot 
imply such a behavior reliably. Loop quantum cosmology, on the other hand, has a dis- 
tinguished ordering and thus predicts the existence of a critical frequency for the system 
studied here. 

Before interpreting the physical significance of uj c we note that the possibility of a non- 
oscillating wave function for uj 2 < uj 2 is not in contradiction with semiclassical behavior 
at large volume. The standard semiclassical analysis requires oscillating (WKB) solutions 
which in our case can be derived from (^) after replacing c with ^ndS/dp and are given 
by N WKB (p) = exp(iS/h) = exp(iQ, WKB hgp) with VL WKB = ^y/3Kuj/lp. For uj 2 > uj 2 the 
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frequency Q in fact reduces to Qwkb- The non-oscillating solutions for u 2 < u 2 cannot 
be interpreted in this way, but this should not be expected since in the classical limit 
uj c = ^kH 2 vanishes and so any non-zero frequency will be larger than the critical one. 



4.2 A critical energy 

We have seen that loop quantum cosmology predicts the presence of a critical frequency for 
a dynamical evolution in the particular model studied here. While this prediction is reliable 
(in contrast to standard quantum cosmology where it may or may not exist, depending 
on the factor ordering of the constraint), we will see that the precise value depends on 
quantization ambiguities of a different kind. 

For an interpretation of uo c we first compute the associated (1/a-dependent) eigenvalues 
of the matter Hamiltonian, 

E(a) = \\a\-*uj 2 . (27) 

In standard quantum cosmology, this expression is unbounded from above and does not 
have a distinguished value which would be suitable for an interpretation. In loop quantum 
cosmology, on the other hand, we do have — for a given u — an upper bound for E which 
can serve as a natural value for an interpretation. Taking for l/\a\ the maximal eigenvalue 
of the inverse scale factor operator, occurring at n = 2, we get a relation between the 
maximal field energy concentrated in Vi/ 2 = |(7^>) 5 an d u 2 : 

oil 



2 

— - UJ 



By inserting uj 2 into ( p5|) the distinction between oscillating and non-oscillating behavior 
can be expressed in terms of this initial energy, 



n Sv/vr 3 lE max 256 _3 

fi = ^r 4 V^~^ 7 2 ' (28) 

with Ep denoting the Planck energy 8nh/lp = Ip/G. This expresses Q in terms of the 
ratio E max /E P , whose critical value ff^T" 5 « 1.006 is close to one for 7 equal to one, but 
has the larger value ~ 21.5 for the value 7 = log2/7rv / 3 ~ 0.13 which has been computed 
by comparing the black hole entropy resulting from a counting within quantum geometry 
with the semiclassical Bekenstein-Hawking result [IS, l!|. Only if E max /Ep is larger than 



the critical value can a dynamical evolution of the classical Friedmann model set in. 

This observation implies that we need a maximal energy in a Planck volume which 
exceeds the Planck energy, being apparently in conflict with the heuristic but widespread 
expectation that there can be at most an energy amount of Ep in a Planck volume. (Often, 
the holographic principle |[2(J is used to arrive at this conclusion. In the present context, 



there is no direct contradiction of results obtained within the formalism used here.) Usually, 
inflation is invoked to explain how the huge amount of energy in the present universe can 
have originated from a Planck scale universe. Here we can see that the relation between the 
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maximal energy and inflation is even more intimate in loop quantum cosmology: One can 
reduce the critical value for E max /Ep by using a different quantization for the inverse scale 
factor instead of (|8]) parameterized by an ambiguity parameter j (a non-zero half-integer, 
see Hl7i) such that 



H7» = 144(j(j + l)(2j + l))- 2 ( 7 ^)^ fc v /Fi (|n+2fe| _ 1) ) . ,20! 

As derived in [fT/jl , the maximal value of H~ 3 is then attained if n = 2j and approximately 
given by V^ 1 (the approximation gets better for large j). This gives a maximal energy 

and 



f2 ~ V47r3 4 (7j)' 



-^p 167r(7j)2 

The critical value ^|(7j) ^ ~ 0.74j~l for E maiX / Ep is then suppressed by j~t and much 
smaller than one for large j (note that this critical value for j = ~ does not coincide with 
the one obtained above because the approximation for the maximum of (|29| ) is bad for very 
small j). Moreover, the maximal energy is not obtained in the initial Planck volume but 
in a volume of size Vj ~ (|j\/7^p) 3 - As observed in [21], using a quantization with a large 



value of j also leads to a prolonged phase of inflation in the very early universe, displaying 
the relation between a maximal energy below the Planck energy and inflation. 

A large value of j modifies the small volume behavior of the wave function considerably 
(Fig. |3|). Its oscillation length and amplitude first decrease with increasing volume, which 
is characteristic of an accelerating universe. When the maximal value for \a\j 3 is reached, 
standard behavior fl26| ) sets in with increasing oscillation length toward larger volume. 



5 Comparison with continuum initial conditions 

In loop quantum cosmology the initial conditions for the gravitational part of the wave 
function are derived from the evolution equation and thus fixed [[H|. While sq, the value 
of the wave function at the classical singularity, drops out of the evolution equation and 
so remains unspecified, the lowest values of n show that the wave function approaches 
the value zero for small n (see Figs, p] and |]). Interpreted as an initial value, this is 
reminiscent of De Witt's proposal that the continuum wave function ip(a) should vanish for 
a = It is well known that DeWitt's initial condition can be satisfied non-trivially 
only in special systems (e.g., for quantum de Sitter space); generically it would imply that 
the wave function vanishes identically. 

An example for a system where DeWitt's condition does not work in general is the one 
discussed in the present paper: (BB) is a continuous wave function which does not have a 
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800 



Figure 3: A wave function with an extended inflationary phase (j = 200) and standard 
behavior at large volume for n > 2j = 400 {uj 2 = 500u 2 ). 

well-defined limit for a — > if uj 2 > uj 2 because it oscillates with constant amplitude and 
diverging frequency close to the classical singularity (for uj 2 < uj 2 De Witt's initial condition 
can be defined and coincides with the result of the dynamical initial condition; see Fig. 0). 



Another example is P2j where a wave function of the same kind as Q2ED was obtained 
for stiff matter and rejected because it does not satisfy a DeWitt boundary condition in 
the sense that either the function or its derivative should be equal to zero at the classical 
singularity. The same is true for the continuum approximation of our wave function, but 
in this case the problem is cured by discreteness near the singularity. 

The oscillations are suppressed once their oscillation length becomes smaller than the 
Planck length, and consequently their amplitude approaches zero. (If there is prolonged 
inflation with a large value of j, as discussed above, the oscillation length may never reach 
the Planck length where this mechanism would be necessary. For a fixed j, however, a 
regime with small oscillation length will always be present if uj is large enough.) This can 
be seen by following the analysis of [1~0||: If we define 



P(n) := iK 7 $i^(n)(V n/2 - V^y 1 
which at large volume is approximately P ~ ^nH^/a, the discrete evolution equation (^D| 
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is 

tn+4-(2-l 2 P(n))t n + t n _ A = 0. 

In the neighborhood of a fixed n we can assume that P(n) is a constant and solve the 
difference equation with an ansatz t n = e me which leads to 

e 4id - (2 - 7 2 P) + e~ m = 2 cos 46 - 2 + 7 2 P = . 



X 

x 




/x 



x 

5 10 15 20 

m 

Figure 4: The small- m behavior of the 0- independent part s^ m [<\))e~ lu) ^l h of the discrete 
wave function (x) compared to a continuous solution to the Wheeler-DeWitt equation 
(dashed line) for uj 2 = 500u 2 . Whereas the continuous solution displays an infinite number 
of oscillations, the discrete solution decreases toward zero once the oscillation length is 
comparable to the Planck scale. 

If P is small, the solution 6 = ^arccos(l — 7 2 P/2) is real and small, implying an 
oscillating wave function t n = e' m6 with long oscillation length 2n/9 ks ir/'jy/P as observed 
at large volume. (This conclusion is only valid if the resulting oscillation length is smaller 
than the length scale on which P(n) changes substantially. In particular, it does not apply 
to subcritical solutions for which the predicted oscillation length would be too large.) As 
P becomes larger, which will happen for decreasing a, will first increase up to # max = 7r/4 
(for P = 2 / ~f~ 2 ) and then become imaginary. The result is that the wave function will 
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first oscillate more and more rapidly, as expected from the continuum approximation, but 
then enter a branch with exponential behavior which does not appear in the continuum 
formulation. Of the two independent solutions — exponentially increasing or decreasing 
- only the decreasing one is allowed by the dynamical initial conditions (Fig. This is 
the mechanism, essentially relying on the discreteness, which allows to generalize De Witt's 
initial condition to systems where the continuum version does not work. (Note that a 
similar behavior with exponential solutions can occur at large volume in the presence of 
a positive cosmological constant. In such a case, it would be caused by large volume 
rather than large curvature and therefore signals an infrared problem. In the small volume 
regime, however, a modified behavior is perfectly admissible and welcome since the classical 
description is expected to break down if curvatures become large.) 

The diverging oscillation number close to a = renders inapplicable not only De Witt's 
proposal but any condition which requires properties of the continuous wave function at 
a = 0, which includes the "no-boundary" and the "tunneling" proposals. (Usually, this 
problem is avoided by setting u = or an analogous condition by hand, see e.g. p3| . 
Note that this is a very strong assumption since even a tiny uj would eventually yield a 
large kinetic term due to the inverse volume.) As already discussed, of the proposed initial 
conditions only the discrete formulation of loop quantum cosmology can deal with the 
problem of wild oscillations close to the classical singularity. 



6 Conclusions 

In this paper we considered a free massless scalar field coupled to loop quantum cosmology 
as a model for implications of the kinetic term in a matter Hamiltonian. This term diverges 
classically at zero volume and is also problematic in standard quantum cosmology where it 
leads to infinitely many oscillations of the wave function close to the classical singularity. It 
is caused by the inverse volume in the kinetic term and is not sensitive to the special form 
of matter or its quantization. Therefore, we restricted our attention to the free massless 
scalar field quantized in a standard way as usual in the Wheeler-DeWitt approach. A 
possible mass or potential term, which is proportional to the volume, would not change 
the qualitative behavior at small enough volume where it would be suppressed. Also higher 
spin fields or different quantization techniques applied to the matter field, e.g. inspired by 
methods for full quantum geometry ||15|| , would not be significant as far as the present paper 



is concerned. We also note that we simplified the analysis by using only the Euclidean part 
of the gravitational constraint multiplied with 7~ 2 instead of the full Lorentzian constraint 
in (|9]). In the flat case, both expressions agree classically, and at the quantum level the more 
complicated Lorentzian constraint does not lead to significant changes of the qualitative 
behavior. 

The quantization of the gravitational degrees of freedom, on the other hand, is impor- 
tant since it affects the form of the inverse volume in a quantization of the kinetic term. 
While standard quantum cosmology treats the inverse volume as a multiplication operator 
which does not cure its divergence, loop quantum cosmology leads to a quantization with 
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significant changes at small scales where the discreteness of the volume is important and its 
inverse does not diverge. This has already been seen to imply the absence of cosmological 
singularities ||, dynamical initial conditions for the wave function of a universe [fTOfl , and 
a new origin of inflation ||21|| . The main result of the present paper is that the discrete- 
ness at small scales also leads to a cure of pathologies of the wave function in a standard 
quantization, like the infinite number of oscillations caused by a kinetic matter term. 
Also the issue of initial conditions is further elucidated by this analysis. In [|10]] it 



has been shown that loop quantum cosmology predicts dynamical initial conditions for 
the wave function of a universe which can be derived from the evolution equation and 
need not be imposed by hand. In this derivation the structure of the cured classical 
singularity plays an important role. In their effect on the wave function the dynamical 
initial conditions resemble most closely DeWitt's initial condition that the wave function 
should vanish at the classical singularity. The drawback of DeWitt's approach, namely that 
this initial condition cannot be fulfilled non-trivially for most systems of physical interest, 
has been seen to be eliminated by effects of the discreteness of loop quantum cosmology. 
Therefore, the dynamical initial conditions of loop quantum cosmology present a meaningful 
generalization of DeWitt's initial conditions. 

We also emphasize that loop quantum cosmology leads to more reliable results than 
standard quantum cosmology because some quantization ambiguities like the factor order- 
ing of the gravitational part of the constraint are fixed. Nevertheless there are ambiguities 
of a different kind which can affect physical consequences and therefore in principle lead 
to observable effects. This has played a role in the discussion of a critical matter energy 
necessary for an oscillating wave function. Its value depends on the quantization ambi- 
guities, but its presence can be concluded from loop quantum cosmology, unlike standard 
quantum cosmology where it may or may not occur depending on the ordering. 
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